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Abstract 

Classical partial oxidation processes often suffer from low selectivities. Promising alternatives are electrochemical 
processes where the oxidation takes place at a packed-bed anode while an oxygen-consuming gas-liquid membrane 
is used as cathode. As a basis for the reliable design of such a process, the performance of oxygen-consuming gas- 
diffusion electrodes (GDE) is investigated experimentally and is analysed based on a rigorous model accounting for 
the reaction microkinetics and all relevant mass and charge transport phenomena. The results indicate that oxygen 
is transported in the gas-filled pores by Knudsen diffusion and that the cathodic oxygen reduction follows a parallel 
reaction scheme forming hydrogen peroxide at the carbon black support and water at the applied platinum catalyst 
particles. 
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mass transfer coefficient 




n 

1 

number of transferred electrons 

ai 0 
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exchange current density with 

P 

Pa 

pressure in gas phase 



respect to electrode volume 

r 

m 

length coordinate through pore 
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dimensionless parameter groups 
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Biot mass number (Table 1) 
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universal gas constant, =8.314 
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mol m -3 
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J mol _1 K“‘ 
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thickness of electrode 
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T 

K 

temperature 
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ionic mobility 
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diffusion coefficient 
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electrode potential vs. standard 
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effective ratio of reactants with 
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charge number 
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1 

charge transfer coefficient 

E 
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enhancement factor of 

p 
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Prater number of electrochemical 



electrochemical reaction 



reaction (Table 1) 
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porosity (void fraction) of 

F 
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Faraday constant, 
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K 
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ionic conductivity of electrolyte 

H 

Pa 

Henry coefficient 

<t> 

V 

electric potential 

Ha 

1 

Hatta number (Table 1) 

V 

1 

stoichiometric coefficient in reaction 

i 

A m“ 2 

current density 


V 

overpotential 
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p 1 electrolyte resistance Lower indices 

a (Q m) -1 ionic conductivity of solid phase 




a 

anodic 



c 

cathodic 



counterion 

counterion in electrolyte 

Upper indices 

eff 

effective 



H + 

protons 

dir 

direct reaction path of oxygen reduction 

i 

of the instantaneous reaction 

e 

in equilibrium 

j 

component index 

FL 

in flooded layer of the electrode (Figure 2) 

i = 0 

at zero current (open circuit) 

G 

in gas film (Figure 2) 

lim 

limiting 

GL 

in gas-filled layer of the electrode 

max 

maximum 


(Figure 2) 

n 

normal to solid surface, towards 

ind 

indirect reaction path of oxygen reduction 


electrolyte 

L 

in liquid film (Figure 2) 

o 2 

oxygen 

S 

in solid phase 

RE 

reference electrode 

* 

dimensionless 

t 

total 


1. Introduction 

In the early 1980s it was suggested that oxygen¬ 
consuming electrodes might be used as counter elec¬ 
trodes in chlorine production, instead of hydrogen- 
producing electrodes [1, 2], This would lead to lower 
energy and safer operation, as hydrogen is highly 
explosive, especially when accidentally mixed with 
chlorine. Other possible applications are as counter 
electrodes in the production of metal solutions (e.g. 
Sn(BF 4 ) 2 ) for electro-plating [3] and in electrochemical 
hydrocarbon partial oxidation processes in the liquid 
phase [4-7]. 

Finally, and maybe most importantly, oxygen cath¬ 
odes play a major role in all types of fuel cells. Also here 
the mechanisms of oxygen conversion (i.e. the electro¬ 
catalysis) and the role of mass transport of oxygen and 
the reaction products (e.g. water) needs to be fully 
understood. Such knowledge is essential for further 
optimisation of fuel cell performance. In this field a large 
number of recent publications are available. For a good 
overview see e.g. [8], 

As a basis for the reliable design of such processes, the 
performance of the oxygen-reducing GDE is investi¬ 
gated experimentally and analysed on the basis of a 
rigorous model which accounts for the reaction micr¬ 
okinetics as well as for all relevant mass and charge 
transport phenomena. 

Figure 1 shows a schematic of such a GDE: Oxygen is 
transported to the electrode surface through the outer 
gas film (step 1), diffuses through a pore (step 2) to the 
gas-liquid interface and is then absorbed in the liquid 
phase (step 3). The absorbed oxygen diffuses through 
the liquid phase (step 4) to the reaction zone (step 5). 
Simultaneously protons are transported from the liquid 
bulk through the outer liquid film (step 7) and the 
flooded part of the pore (step 6). In the reaction zone 
(step 5), oxygen, protons and electrons (step 8) react to 
form water. As water is the main component of the 


liquid phase, its transport out of the pore has little 
impact on the other transport processes. Therefore, it is 
neglected in the further analysis. 

The GDEs used for the experiments consisted of 
platinum catalyst supported on carbon black (total 
platinum loading 1.5 mg cm -2 ), carbon black without 
catalyst loading and a polymeric binder (PTFE). A 
detailed description of the production process as well as 
information on the physical properties of the GDEs can 
be found in [9, 10]. 


2. Process analysis 


In the following, a mathematical model for the 
simulation of polarisation curves of GDEs will be 
derived, based on an approach previously developed 
for first-order reactions [11]. In order to obtain a 
compact model formulation, dimensionless variables 
and parameter groups are introduced, marked with 
the upper index They are listed and explained in 
Table 1. 

The electrochemical oxygen absorption reaches a 
maximum value of the current density, 4ff, max 
[A m -2 ], when only the oxygen diffusion in the gas filled 
part of the pores (the gas layer, step 2 in Figure 1) is the 
rate limiting step: 


k-fT.n 


= nF 


:\’G[- n GL G 

iV 02,max _ PO 2 

vq 2 ~ n cP L RT 


( 1 ) 


Here n [-] is the number of transferred electrons per 
net reaction, F— 96485 C mol -1 is Faraday’s constant, 
7Vo 2 L ma x [ m °l m 2 s_1 ] the maximum oxygen molar flux 
density, [m 2 s -1 ] is the effective oxygen diffusion 

coefficient in the gas layer (GL), d° L [m] is its thickness, 
Pg 2 [Pa] is the oxygen partial pressure in the gas film (G), 
R = 8.314 J mol -1 K -1 is the universal gas constant 
and T[K\ is the gas temperature. 
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Fig. 1. Structure of gas diffusion electrode (GDE) and transport phenomena occurring for oxygen absorption. 


Table 1. Definition and interpretation of dimensionless parameter groups 


Biot mass number for gas phase 

Biot mass number for liquid phase 

Thickness of flooded layer (FL) of electrode 

Ratio of mass transport resistances (GL/FL) 

Effective reactant ratio with respect to bulk phases 

Enhancement factor of electrochemical reaction 

Enhancement factor of instantaneous reaction 

Flatta number 

Effective current density 

Prater number of electrochemical reaction 

Electric resistance of electrolyte 


Definition 


BE 3, 


Bib 3 


k % 2 d 

JgT— 

eff, 0 2 
kk d 


L _ k o 2 


d* = 


= 

d 


= B So 

_ c Q 2 (Pq 2 )/ V 02 

e 4+/%+ 


E = - 


E,= 1 + 


Ha = d 


- ss Her,Ha = d*Haexp(— 0.5a c f/*) 

~FL , 
v Q2 £) eff,H+ c H + 

v h+ u So 2 ‘o>8,) 

aio/(nF) 

171 c Le (p G ) 


P = 


left, max 

_" F2 °5o 2 ‘^(t > g 2 ) 


aiJtT 


, _ 0£ O2 pS 2 d " F ~ 

p — >’kv'(itr 


Physical meaning 


mass transfer in G 
pore diffusion velocity in GL 

mass transfer in L 
pore diffusion velocity in FL 

total electrode thickness 

pore diffusion velocity in FL 
pore diffusion velocity in GL 

0 2 concentration at liquid phase 
FLconcentration in L 

Q 2 absorption with reaction 
0 2 absorption without reaction 

Fl + pore diffusion velocity in FL 
0 2 pore diffusion velocity in FL 

reaction rate with respect to L 
pore diffusion velocity in FL 

effective current density 
maximum effective current density 

pore diffusion velocity in FL 
electron transport in S 

pore diffusion velocity inGL 
charge transport in L 


Total overpotential 


— rt(U U}=o) 


This maximum effective current density is used to 
formulate the dimensionless current densities according 
to Table 1. 

2.1. Model for electrochemical gas absorption in gas 
diffusion electrodes 

2.1.1. Mass transport 

Figure 2 shows schematically the qualitative reactant 
profiles within a pore of the electrode for low current 
densities (i.e. slow reaction). As the main transport 
resistances are expected within the flooded layer of the 
pore, the concentration profiles within the outer gas and 


liquid film, respectively, as well as within the gas layer 
within the electrode are assumed to be linear. The 
following expressions describe the molar flux densities 
Nj [mol m -2 s -1 ] in these layers, assuming Fickian 
diffusion: 

Liquid film (L): Nf s = kf ( C/ L - c} s ); j = 0 2 , H+ 

( 2 ) 

p G - p GS 

Gas film (G): Nf s = k f 1 ^ ■ j= 0 2 (3) 
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Fig. 2. Qualitative oxygen and proton concentration profiles in 
GDE. 


Gas layer (GL): Nf L 


RT ’ 


j = o 2 

(4) 


In these equations, Cj [mol m -3 ] are molar concentra¬ 
tions in the liquid phase, pj [Pa] are partial pressures in 
the gas phase, and kj [m s -1 ] are mass transfer coeffi¬ 
cients. 

Within the gas layer, mass transport is not only driven 
by molecular diffusion, but also by Knudsen diffusion. 
The superposition of both mechanisms can be described 
by an effective diffusion coefficient D^ f h 


1 


D 


GL 

effy 


1 -pf/p° 1 
n a + D K 

V e((j U eff.j 


(5) 


calculated from the effective diffusion coefficients in the 
(free) gas phase, and the effective Knudsen 

diffusion coefficient -D£j Ty . 

The gaseous reactant (oxygen) dissolves in the liquid 
phase at the gas-liquid interface. This process is 
assumed to be in equilibrium, and the concentration 
on the liquid phase side, c Fe , is calculated from the 
partial pressure on the gas phase side, pj, according to 
Henry’s law: 

pf = Hj C y/c\ (6) 

where Hj [Pa] is the Henry coefficient and c\ is the total 
concentration in the liquid phase. 

Within the flooded layer, mass transport can take 
place due to diffusion and electro-migration. The latter 
is driven by the gradient of the electric potential, cj> FL 
[V], and only affects charged species: 


N] 


rFL = D 


ef 


dr 


Here Zj [-] is the charge number of species /', . 

[mol nr J -1 s' 1 ] is its ionic mobility and r [m] is the 
spatial coordinate (see Figure 2). 

Equation (7) is only valid for dilute electrolytes. To 
account for the porosity e [-] of the electrode, the 
following approximation for the molecular transport 
coefficients were used (as proposed e.g. by Newman 
[ 12 ]): 

n G „1.5 n G n FL „1.5 r)L FL _ .1.5, .L ( o\ 

U eff,/ ~ s U j ■ U ef(j U j ' M eff J ~ s u j l 8 ) 

For the model, water is not accounted for as it is the 
excess component in the liquid phase and should 
therefore have no significant influence on the overall 
process. 

Finally, in the flooded layer mass balances for oxygen 
and protons, as well as the charge balance, are formu¬ 
lated in terms of the dimensionless numbers given in 
Table 1: 

°xyg en: -j^r= -(d*Ha) 2 i* (9) 


Protons: 


d 2g H + _ Ad* Ha) 2 
dr* 2 Ej - 1 " 


( 10 ) 


The factor 4 in Equation (10) results from the 
stoichiometry of the electrochemical oxygen reduction: 

0 2 (g) + 4e~ + 4H + —> 2H 2 0 (11) 


Ej [-] stands for the enhancement factor of the 
instantaneous reaction (i.e. maximum current density, 
reaction takes place in a planar reaction zone within the 
flooded layer). 

The dimensionless boundary conditions for the oxy¬ 
gen balance are: 


dc7 


dr* 


1-c* (r* = 0) 


D* 

Bi °d> 




( 12 ) 




L o, 


(r* = 1) - cj 


,*L 

o 2 


1 


(13) 


Pore diffusion in the gas layer is described by D*, 
the mass transport resistances in the outer liquid and 
gas film are accounted for by the respective Biot 
numbers, Bi^ and Bi F . In the same manner, boundary 
conditions can be formulated for the proton mass 
balance: 


dr* 

af H+ 

dr* 


r*=o 


= 0 


(14) 


dr* 

aC H + 

dr* 


BinX 

e{E,~ 1) 


( 4 ^ 


1 )- 1 ) 


(7) 


r *=1 


(15) 
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with e being the effective ratio of both reactants related 
to the bulk phase conditions, taking into account the 
respective concentrations and stoichiometric coeffi¬ 
cients: 


_<&(?&)/(" 1) 

c h+/(—4) 


(16) 


2.1.2. Charge transport 

Charge transport takes place in the solid matrix of the 
GDE (electron transport) as well as in the electrolyte 
(ion transport). Within the solid matrix, charge trans¬ 
port is assumed to follow Ohm’s law 


r = —a 


eff 


d</r 

dr 


(17) 


The effective conductivity of the porous matrix is 
approximated by 

«4=(l-e) 1 V. (18) 


Charge transport within the electrolyte can be 
described by Faraday’s law accounting for the fluxes 
of all charged species: 


N 

i FL = zjNj L ; / = H + , counterion (19) 

7=i 


The charge balance for the solid matrix can be 
formulated as 


— = -ai„ (20) 

where a [m 2 m“ 3 ] is the volume-related inner electrode 
surface, and i n [A m~ 2 ] is the local current density 
normal to the inner solid surface. 

For the electrolyte in the flooded electrode layer, 
electroneutrality is assumed: 

X>/f = 0- (21) 

7=1 

This assumption is only strictly valid for the core of 
the electrolyte phase, not for the electrochemical double 
layer. But it is justified, as long as the total volume of the 
double layer is small compared to the total volume of 
the pore liquid. However, for highly porous systems, 
significant deviations from electroneutrality can occur. 
The same is true for highly dilute electrolytes, where the 
double layer thickness can reach values of more than 
10 nm. 

For a binary electrolyte, accounting for the stoichi¬ 
ometric factors s H + and 5 counte rion of the two ionic 
species within the electrolyte molecules, the molar 
electrolyte concentration is given by 


This equation couples the concentration profiles of the 
two ionic species in the electrolyte and therefore the 
concentration of the non-reactive of the two species 
(counterion) can be eliminated. Its molar flux density is 
constant throughout the pore and equal to the value at 
the gas-liquid interface. For the molar flux densities of 
the two ionic species, one consequently gets the follow¬ 
ing expression: 


«f L = 


FL 


J eii,2“y 

j = H + , counterion 


dr 


(23) 


From this, the profile of the electrical potential 
(concentration overpotential) can be derived as a 
function of the concentration profile c FL (r): 

d<j> _ "®eff,counterion d ^ ^f-j 1 ( 24 ) 

dr “counterion Meftcounterion^ dr 


Using this equation, it is possible to calculate the 
molar flux of the protons solely from their concentration 
gradient: 


N 


FL 

H+ 


= -D 


FL 

eff,H+ 



with the effective diffusion coefficient 


(25) 


D 


FL 

eff.H+ 


^counterion ^eff counterion ^eff.H 1 “ ^ ^ W eff,H + ^eff,counterion 


rFL 


FL 


~>FL 


“counterion w dT, C ounterion 


(26) 


which characterises the whole electrolyte. 

Finally, inserting Ohm’s law, Equation (17), into the 
charge balance of the solid matrix. Equation (20), one 
obtains in dimensionless notation 


d 2 A<j> s * 

dr* 2 


-fid* 2 Ha 2 ? 


(27) 


The parameter (3 describes the influence of the electric 
resistance of the solid matrix in relation to the pore 
diffusion resistance of oxygen and is an analogue to the 
Prater number known from heterogenous catalysis. 
Therefore, it is justified to call (3 the Prater number of 
the electrochemical reaction (see Table 1). 

The dimensionless potential profile of the pore liquid 
can be obtained from Equation (24) as 


dA <f> v 
dr* 


= —D 


* 

counterion 


d In 
d r* 


(28) 


The diffusional and migrational properties of the 
counterion are lumped in the parameter £>*^ terion . For 
very dilute electrolytes the ion mobility of the counter¬ 
ion can be calculated as 


c FL r FL 

^FL _ H + _ ^counterion 

^H + ^"counterion 


( 22 ) 


u t 


FL 

eff,counterion 


d fl 

eff,counterion 

RT 


d : 


^counterion 


( 29 ) 
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In the following, this relation will also be used for 
electrolyte solutions with medium concentrations as a 
fair approximation. 

For the boundary conditions for the charge trans¬ 
port in the solid matrix, only the ohmic resistance in 
the matrix has to be accounted for, as long as the 
resistance in the current collectors (see Figure 1) can 
be neglected. Charge transfer reactions only occur at 
the pore walls within the GDE, but not on the GDE’s 
outer surface: 


dA(f> s * 

dr* 


r*= o 


1 -d' 


■A 0 s *(r* 


dA<(> s * 

dr* 


= 0 


;■*=1 


0 ); 


(30,31) 


Finally, for the charge balance of the pore liquid an 
initial condition is needed. Here it has to be taken into 
consideration that the potential difference Ac|) FL is given 
with respect to a reference electrode RE (see Figure 2), 
which defines the potential within the liquid bulk. 
Between there and the edge of the boundary layer 
(distance d RE ), due to the electrolyte resistance, a 
potential difference proportional to the electric current 
density occurs (7-R-drop) 

(<£ L ) L - (<^ L ) RE = —(32) 

The total potential difference between electrode sur¬ 
face and reference electrode is the sum of the i-R-drop, 
Equation (32), and the concentration overpotential in 
the liquid film. The latter is obtained from integration of 
Equation (28). Putting all this together, one gets in terms 
of the dimensionless numbers: 


A(/> l * (r* = 1) = In (ch+(r* = 

VCh 

(33) 

with the dimensionless electrolyte resistance p~. The sign 
of the stoichiometric factor v 02 determines whether an 
anodic or cathodic current is flowing. 


2.1.3. Reaction microkinetics 

For the local current density of the electrochemical 
reaction, a Butler-Volmer micro kinetic rate expression 
is applied. In dimensionless notation, it is given by: 

i* n = exp{ot a (rf - A</> s * - A(p L *)} c*"“ J 

J=l n (34) 

- exp{—a c (//* - A(j) S * - A0 l *) } cf J 

7=1 

where a a and ot c [-] are the charge transfer coefficients of 
the anodic and the cathodic reactions, respectively, and 
r| [-] is the dimensionless reaction overpotential. 

The combination of Equations (34), (9) and (10) yields 
an expression for the effective current density: 


* *1 -d* 1 _ 1 -d* 

,e ff- ct RG + R GL + R FL- 1 +( y _#)+£_■ 

(35) 

The three terms in the denominator, R°, R ( l and 
R fl , are the mass transport resistances in the gas film, 
the gas layer and the flooded layer, respectively, as well 
known from classical absorption theory. 

To get analytical solutions of the set of differential 
equations, the concentration profiles in the liquid phase 
are approximated by piecewise linear profiles as pro¬ 
posed by van Krevelen and Hoftijzer [13, 14]. This leads 
to the following expression for the enhancement factor 


tanh Hao 2 


with 


(36) 


Ha = d* Ha exp y (r* = 0) 

= ct Ha exp y rf^j Q. _f c H+( r * = 

(37) 

The latter is the modified Hatta number for electro¬ 
chemical oxygen absorption, formulated with respect to 
the proton concentration on the gas side of the GDE, 
Cy + (r* = 0), first line, and with respect to the proton 
concentration on the liquid side, c^ + (r* = 1), second 
line (positions: see Figure 2). 

In Figure 3, simulated polarisation curves are plotted 
according to Equation (35) under variation of the 
enhancement factor for the instantaneous reaction E f . 
All parameters used are given in Figure 3. Three distinct 
operating regimes are evident - slow, fast and instan¬ 
taneous reaction - for which qualitative concentration 
profiles are also depicted. For low overpotentials (slow 
reaction), the oxygen reduction kinetics are slow and 
limit the overall current density. A mass transport 
limitation can be neglected and the reaction takes place 
throughout the flooded layer (oxygen concentration 
profile is decreasing throughout the flooded layer). For 
higher overpotentials, the reaction rate increases (fast 
reaction) and the overall process is now increasingly 
affected by internal mass transport resistances. There¬ 
fore in this regime the polarisation curves show a clear 
dependence on the E-factor. The oxygen concentration 
profile lies much lower than in the previous regime and 
the oxygen concentration reaches zero within the 
flooded layer, i.e. the absorbed oxygen is fully con¬ 
sumed. Obviously, the reaction takes place only in a 
zone near the gas-liquid interface. This zone is the 
smaller, the higher the overpotential (and therefore the 
overall reaction rate) is. A further increase in overpo¬ 
tential leads to an instantaneous reaction of oxygen and 
protons. The reaction zone is compacted to a reaction 
plane, the process is limited by mass transport of the 
reactants oxygen and protons only. The different mass 
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Fig. 3. Simulated polarisation curves according to Equation (35) and qualitative concentration profiles for three different operating regimes. 


transport mechanisms determine the location of this 
reaction plane within the flooded layer. The correspond¬ 
ing current density reaches a limiting value which only 
depends on the enhancement factor of the instantaneous 
reaction E h The higher this enhancement factor, the 
higher the effective limiting current density. For an 
infinite enhancement factor (E t —> °°), an asymptotic 
polarisation curve is reached, which corresponds to a 
pseudo first-order reaction with respect to oxygen. 

2.2. Experimental results 

To validate the model and to identify its parameters, 
potentiostatic polarisation experiments were performed 
using the GDEs described in Section 2. The experiments 


were carried out in a specially designed cyclone flow cell. 
This provides for a well-defined flow field on both sides 
of the GDE (comparable to a rotating disk arrange¬ 
ment) and allows the independent adjustment of the 
outer gas and liquid flow conditions. A detailed descrip¬ 
tion of this cell is given in refs. [15] and [10]. 

In the experiments, pure oxygen was used. The oxygen 
flow rate, the total gas pressure and the electrolyte 
concentration were varied as most important process 
parameters. The results show that, of these, the electro¬ 
lyte concentration is most significant. Figure 4 presents 
some experimental polarisation curves for different 
electrolyte concentrations. From the model, one would 
qualitatively expect curves like the dashed one, if a single 
electrochemical reaction is involved. However, the 
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Fig. 4. Experimental polarisation curves for different electrolyte concentrations and qualitative comparison with model Equation (35) (dotted 
curve). 
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experimental curves, show a significantly different 
behaviour. Going from right to left (i.e. from low to 
high overpotentials), first the current densities increase, 
but seem to reach an intermediate limiting value, before 
finally increasing further. This indicates a reaction 
mechanism which is more complicated than anticipated 
in the model at this stage. 

2.3. Analysis of oxygen absorption 

According to the literature (e.g. [16]), oxygen reduction 
in an acidic environment can follow two possible 
reaction paths: A direct one (upper index “ dir ”) 

0 2 (g)+4e-+4H + —>2H 2 0 if = 1.229 V (38) 

and an indirect one (upper index “ind”) 

O 2 (g) + 2e“ + 2H+ -f H 2 0 2 if = 0.623 V 

H 2 0 2 + 2e“ + 2H+ -> 2H 2 0 U 9 = 1.834 V 

(39,40) 

being a consecutive reaction with hydrogen peroxide as 
intermediate. The direct path is catalysed by platinum, 
palladium, silver and iron-N 4 -chelates and has a stan¬ 
dard cell voltage of = 1.229 V. The first step of the 
indirect path can be observed on carbon surfaces like 
graphite and carbon black without any metal catalyst, 
and it has a fairly low standard potential of only 
t/l“g = 0.623 V. The hydrogen peroxide decomposition 
in the second step is catalysed by nickel and cobalt and 
has a very high standard potential of 1.834 V. 

Based on these phenomena, the following assump¬ 
tions have been made for the modification of the 
mathematical model: 

(a) On the platinum catalyst, oxygen reacts to water 
using the direct path (4-electron-transfer). 

(b) On the carbon black surface, oxygen reacts to 
hydrogen peroxide (2-electron-transfer). 

(c) Hydrogen peroxide decomposition is slow com¬ 
pared to its formation and can be neglected [17]. 

(d) Far away from the equilibrium potential of the 
indirect oxygen reduction (lf = 0.623 V), the ano¬ 
dic backward reaction can be neglected. 

(e) Both reactions microkinetics are described by a Ta- 
fel approach with the reaction orders being identi¬ 
cal to the stoichiometric coefficients: 


Direct path: iff = -Cq 2 (c h+ *) 4 



Indirect path: 4ff = ~ c (>>( <: 'h +) 2 exp{-a[ IK V} 

(42) 

The overpotential q* is formulated with respect to 
the indirect reaction path (see Table 1), therefore in 
Equation (41) the term — (7™$) appears to 

account for the overpotential of the direct path. 


(f) The direct oxygen reduction is slower than the indi¬ 
rect reduction, despite the presence of the platinum 
catalyst (some orders of magnitude higher carbon 
surface than platinum surface). 

Based on these assumptions, the dimensionless effec¬ 
tive current densities of both reaction paths are derived. 


2.3.1. Direct oxygen reduction path (upper index “dir”) 
According to Equation (35) the cathodic current density 
of the direct reduction path can be formulated as 


•dir* _ 

*eff — 


l -<f 


(1 “ d *) + WE^ 


D*- -— £ dir . 

d* 


(43) 


As in our experiments pure oxygen is used, no 
transport limitations occur in the gas phase outside the 
electrode and consequently the Biot number Bi(^ is 
infinite. Therefore the term in Equation (43) connected 
with this Biot number vanishes. According to assump¬ 
tion (f) (see list in previuos section) that the direct path 
only has a moderate reaction rate, the enhancement 
factor according to Equations (36) and (37) can be 
approximated by 

Z7 1,r w Ha du = (fHa dir 



(44) 


As already mentioned in the previous section, Equa¬ 
tions (41) and (42), the overpotential q is formulated 
with respect to the indirect reaction path, therefore in 
Equation (44) the term A, ( U d f () — T(=o) * s necessary to 
obtain the overpotential of the direct path. 


2.3.2. Indirect oxygen reduction path (upper index “incT ’) 
Following the same argument used for the direct 
reaction path, the cathodic current density of the 
indirect reduction path is 


•ind* 

*eff 


1 -d* 

(1 ~ct) + 1 ^ rA 

1 -d* 




,* 1— 7> nd * 

a* 1 eff 

D* £j nd -l 


1 


( 

a ind \ 

D*Ha md exp I 



(45) 


with E‘ nd being the enhancement factor of the instanta¬ 
neous reaction with respect to the concentrations in the 
liquid bulk and the gas bulk (see explanation in 
Table 1). In the denominator of Equation (45) three 
resistance terms are summed up (from left to right): 

• the diffusion resistance of oxygen in the gas-filled 
pores of the electrode, 

• the diffusion resistance of the protons in the liquid- 
filled pores of the electrode and 

• the combination of the resistances of the electro¬ 
chemical reaction and the diffusion of oxygen in the 
liquid phase. 

For high overpotentials q the latter resistance 
vanishes and the current density reaches one of two 



possible limiting values, depending on the enhancement 
factor: 
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.•ind 

•ind* eff.lim 

'eff.lim — -ind 

*eff,max 


r)* 1 —d* rind 

d* i 


for Ef<T4=, 
for P nd > 


~D>(\-d>) J 

(46) 


This dependence is shown in Figure 5. For compar¬ 
ison three effective limiting current densities are plotted. 
These were calculated from a lit of the model for the 
indirect path. Equation (45), to the experimental polar¬ 
isation curves given in Figure 4, for low overpotentials. 
The model shows good agreement with the experiments. 


2.3.3. Formulation of the overall current density 

The overall current density is given as the sum of the 

current densities of both reaction paths: 



(47) 



Fig. 5. Effective limiting current density as function of the enhance¬ 
ment factor of the instantaneous reaction, comparison between mod¬ 
el, Equation (46), and values calculated from experimental data. 


with B,=2^ m&x Df\~d*)Hc^ (48) 

and B 2 = ,-' m ax D f \~ d *) Ha md (49) 

In the experiments, the overall current density is 
measured as a function of the electrode potential. 
According to Equation (47) the following five parameter 
groups determine this dependence: 

• the group B x representing the reaction rate of the di¬ 
rect oxygen reduction, 

• the group B 2 representing the reaction rate of the 
indirect oxygen reduction, 

• the limiting current density of the indirect oxygen 
reduction 4 d lim , 

• the cathodic charge transfer coefficient of the direct 
oxygen reduction a dir , 

• the cathodic charge transfer coefficient of the indi¬ 
rect oxygen reduction a : c nd . 

Their values are given in Table 2, determined by non¬ 
linear regression of experimental polarisation data for 
three different proton concentrations in the liquid phase. 
The simulation results show excellent agreement with 
the experimental observations, as can be seen in 
Figure 4 (solid curves). 

The superposition of both reaction paths contribu¬ 
tions to the overall current density is shown in Figure 6. 
For low overpotentials, only the indirect path contrib¬ 
utes to the cell current and finally reaches its limiting 
current density. For higher overpotentials direct oxygen 
reduction also occurs and dominates the process for 
high overpotentials. Within the plotted region, it does 
not reach its limiting current density. The reason why 
the limiting current density of the indirect path lies much 
lower than that of the direct path is most likely the slow 
decomposition of hydrogen peroxide combined with the 
strong adsorption of hydrogen peroxide on carbon 
surfaces. Therefore it can be assumed, that under steady 
conditions most of the carbon surface is covered with 
hydrogen peroxide. 

2.3.4. Oxygen transport in the gas-filled pores 
With the help of the identified parameters (Table 2) and 
Equations (1) and (45), the maximum limiting cur¬ 
rent density of the indirect oxygen reduction can be 


Tabte 2. Enhancement factors E j nd and model parameters for three selected proton concentrations derived from non-linear regression (Confi¬ 
dence intervals valid for a confidence number of 95%), and ratio of exchange current densities calculated from Equation (55) 


Parameter 

= 0.1 mol dm 3 

c[j + = 0.2 mol dm 3 

c[j + = 2.0 mol dm 3 

E ini 

293 

585 

5836 

-B t / A cm -2 

0.93x10“ 6 ±0.11x10“ 6 

1.12x10“ 6 ±0.09x10“ 6 

6.48 x10 -8 ±1.67 x10 -8 

-B 2 / A cm -2 

1.64x10“ 2 ±0.04x10“ 2 

3.07x10“ 2 ±0.05x10“ 2 

4.11x10 -3 ±0.18x10 -3 

-'“lim / A Cm “ 2 

0.0173 ±0.0003 

0.0381 ±0.0005 

0.1016 ± 0.0015 

(Y dlr 

0.48 ±0.01 

0.50 ±0.01 

0.58 ±0.02 

4 nd 

0.76 ±0.02 

0.69±0.02 

0.67 ±0.01 

(mo) di 7(«b) ind 

1.7xl0“ 9 

0.6x 10 _I ° 

1.2xlO“ 10 

(ai 0 ) md /Am 3 

5.7x 10 10 

19.9x 10 10 

3.6x 10 9 
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Fig. 6. Simulated polarisation curve using the final model. Equation 
(47): Superposition of both reaction paths of the oxygen reduction. 


calculated: 


= 0.1016 A cm . From this it is 


I i m d I 
| *eff,max | 

possible to estimate the diffusion coefficient of oxygen in 
the gas filled pores of the electrode (gas layer), assuming 
that the penetration of the liquid phase is small 
(d FL «d) due to the hydrophobicity of the electrode 


matrix: 

rjGL _ D off\(h 

°> ~ el-5 


,-ind 

*eff,max 
~G 


d 


2^p$F/(RT) 


2.77 x Kr 7 m 2 s _1 


(50) 


This value indicates that within the electrode oxygen 
is mainly transported via Knudsen diffusion. The pore 
equivalent diameter for oxygen can therefore be calcu¬ 
lated from the molar mass of oxygen, M Q2 [kg mol -1 ], 
using Knudsen’s equation 

dp.eq = \ £>§ 2 l o, = 1.9 nm (51) 


Obviously, oxygen is transported through the microp¬ 
ores of the electrode. 

Finally, by rearranging Equation (46) one obtains an 
expression for the dimensionless thickness of the flooded 
layer as a function of the enhancement factor of the 
instantaneous reaction E'f d and the effective limiting 
current density of the indirect path 4fr lim : 

D* 

cT= -js- (52) 

7A* _i_ 1 eff.lim 

i~ ^pind ^ind 

i eff,max 


According to Equation (46), the maximum current 
density i") d |im = 4ff max reac hed for enhancement 
factors of the instantaneous reaction of E'fi d =1613 
and higher. Inserting this into Equation (52), the 
thickness of the flooded layer can be calculated: 
d FL D* 

ct=^=- -—r = 0.225 (53) 

d D* + \ 1613 v ' 


with D* 


D°ko 2 PoJ( RT ) 


= 1.8 x lO ^ 4 


(54) 


Therefore, under the given experimental conditions, 
the liquid phase only penetrates less than a quarter of 
the electrode thickness, approximately fulfilling the 
assumption of low penetration. 


2.3.5. Analysis of exchange current densities 
From the values of the parameter groups B\ and B 2 
(Table 2) the ratio of the exchange current densities of 
the direct and the indirect oxygen reduction can be 
calculated as 


(m 0 ) dir _ (Ha dir ) 2 4F _ 1 /5,\ 2 
(az' 0 ) ind ~~ {Ha md ) 2 2F~ 2 \B 2 ) 


The values are also given in Table 2, as well as the 
exchange current densities of the indirect oxygen reduc¬ 
tion, calculated from the parameter group By. 


(m 0 ) ind = 


2FDl (Wl c%M 2 m) 2 


(56) 


It can be seen from these results, that the exchange 
current densities of the direct oxygen reduction to water 
are much smaller than those of the indirect reduction to 
hydrogen peroxide. This has two reasons: on one hand, 
the activity of platinum as catalyst for the direct oxygen 
reduction is low. On the other hand, the platinum 
surface is small compared to that of the carbon black 
surface. 


2.3.6. Analysis of charge transfer coefficients 
Also given in Table 2 are the charge transfer coefficients 
of the direct and the indirect reaction path of the 
electrochemical oxygen reduction, a dir and a[ nd , for three 
different proton concentrations. The values for the 
direct path, a dir , are very close to the ideal theoretical 
value of 0.5 which leads to the conclusion that the 
derived Butler-Volmer-type reaction rate expression. 
Equation (43), describes the process very accurately. The 
charge transfer coefficients of the indirect path, a[. nd , 
have slightly higher values of around 0.7. This leads to 
the conclusion, that the derived rate expression, equa¬ 
tion (45), is also a reasonable description of the 
underlying process, but the approximation is not as 
good as for the direct path. Obviously some of the 
assumptions made for formulating the model with 
respect to the indirect path have to be more closely 
checked. Further experimental results not mentioned in 
this paper show, that there is a slight influence of the 
flow velocity in the liquid bulk. The higher the liquid 
flow velocity, the lower are the current densities. This is 
plausible, considering the influence of peroxide decom¬ 
position, which has been neglected during the derivation 
of the model. The higher the liquid flow velocity, the 
more peroxide is transported away from the electrode 
surface. Consequently, less peroxide is available for the 
decomposition reaction and therefore the contribution 
of this reaction to the overall current is reduced. 
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3. Conclusions 

It can be concluded that oxygen transport in the gas- 
filled part of the GDE pore structure is controlled by 
Knudsen diffusion in the micropores. Due to the 
hydrophobicity of the electrode material, the liquid 
phase only penetrates the macropores with a total 
penetration depth of less than 25% of the electrode 
thickness. 

Furthermore, it was confirmed that the electrochem¬ 
ical oxygen reduction follows a parallel reaction scheme 
forming hydrogen peroxide at the carbon black support 
(indirect path), and water at the platinum catalyst 
particles (direct path), with the indirect path being the 
dominating reaction for low overpotentials. 

Even though hydrogen peroxide decomposition was 
neglected for this first model approach, good agreement 
with the experimental data was achieved. 

With the help of the proposed model for the transport 
and reaction processes, proton concentrations can be 
estimated, for which a maximum current density is 
achieved. This is the key for the estimation of the most 
important electrode parameter: the diffusion coefficient 
of the gaseous reactant in the gas-filled part of the 
electrode. 

Finally, it can be concluded that the derived analytical 
model equations are applicable for the design of 
electrochemical reactors for gas absorption based on 
GDEs. Also the classical film theory of chemical 
absorption is extendable to electrochemical absorption 
processes. 
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